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(57) ABSTRACT 

A method and associated system for use of statistical param- 
eters based on peak amplitudes and/or time interval lengths 
and/or depolarization-repolarization vector angles and/or 
depolarization-repolarization vector lengths for PQRST elec- 
trical signals associated with heart waves, to identify a per- 
son. The statistical parameters, estimated to be at least 192, 
serve as biometric indicia, to authenticate, or to decline to 
authenticate, an asserted identity of a candidate person. 
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Compute at least two of differences A 1(k1), A2(k2), A3(k3) 
and A4(k4) and compare these with the corresponding 
range(s) Rl(kl), R2(k2), R3(k3) and R4(k4), respectively 
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Compute four differences Al(kl), A2(k2), A3(k3) 
and A4(k4) and compare these with the corresponding 
range(s) Rl(kl), R2(k2), R3(k3) and R4(k4), respectively 
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HEART ELECTRICAL ACTIONS AS 
BIOMETRIC INDICIA 

ORIGIN OF THE INVENTION 

5 

This invention was made, in part, by one or more employ- 
ees of the U.S. government and may be made, used and/or 
sold by the U.S. government without payment of royalties or 
other compensation thereon. 

10 

FIELD OF THE INVENTION 

This invention relates to use of a physiological parameter 
as a biometric indicium for authentication of identity of a 
person. 15 

BACKGROUND OF THE INVENTION 

One or more biometric indicia, such as fingerprints, voice- 
prints, retinal scans and facial features, are often proposed to 20 
be used to identify, or to authenticate the asserted identity of, 
a user who seeks access to a given resource. Approximately a 
dozen different biometric indicia have been proposed, but 
implementation methods for some of these approaches have 
not been disclosed. Many of these biometric indicia are asso- 25 
ciated with inherent physiological characteristics of the 
user’ s body. Another set of such indicia relate to what may be 
characterized as neuro -physiological (“N-P”) characteristics 
that partly reflect a learning or behavioral process and do not 
rely exclusively on purely physiological features. Use of one 30 
or more of these N-P characteristics as a biometric indicium 
has received relatively little attention, in part because of the 
perceived difficulty of implementing a procedure to measure 
such a characteristic. An example is a sequence of bioelectric 
signals associated with cycles of the heart. 35 

Cardiac muscle is myogenic and is capable of generating 
an action potential and depolarizing and repolarizing signals 
from within the muscle itself. An intrinsic conduction system 
(ICS), a group of specialized cardiac cells, passes anelectrical 
signal throughout the heart. The ICS includes a sino-atrial 40 
(SA) node, an atrio -ventrical (AV) node, the bundle of His, 
right and left bundle branches, and the Purkinje fibers, as 
illustrated in FIG. 1. These components spread the depolar- 
ization waves from the top (atria) of the heart down through 
the ventricles. The autonomic nervous system modulates the 45 
rhythm, rate and strength of cardiac contraction. When the 
cardiac muscle fibers contract, the volumes within the two 
atrial or two ventricle chambers are reduced and blood pres- 
sure increases. The (smaller) atrial chambers receive blood 
from the veins and pump the blood into the (larger) ventricle 50 
chambers, which pump blood out into the major arteries. The 
heart cycle normally begins in the right atrial chamber, and 
spreads to the left atrial chamber and to the two ventricles. 
The atrial contraction is followed by the ventricular contrac- 
tion in each cycle. 55 

Simultaneous contraction of the large number of fibers in 
the ICS generates an electrical field that can be measured at 
the body surface using an electrocardiograph (ECG). This 
electrical signal includes a sequence of PQRST complexes, 
one of which is schematically illustrated in FIG. 2. The time 60 
interval between two consecutive R signal peaks, referred to 
as an R-R interval, corresponds to a heart pulse, with a rate 
that normally lies in a range of 60-90 beats per minute (bpm). 
The P signal corresponds to atrial depolarization (right side, 
followed by left side) ; the larger QRS complex corresponds to 65 
depolarization of the ventricles and (smaller magnitude) 
repolarization of the atria; and the T signal corresponds to 
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repolarization of the ventricles. A weaker U signal occasion- 
ally appears on the chart, representing remnants of ventricular 
repolarization, but is not shown in FIG. 2. 

According to naming conventions accepted by most cardi- 
ology workers, a time increment with a straight line ampli- 
tude extending between two consecutive signals, for 
example, from the end of a P wave to the beginning of an 
immediately following Q wave, is referred to as a “segment;” 
and a time increment that includes at least one wave, with a 
graph that is at least partly curved, for example, from the 
beginning of a Q wave to the beginning of an S wave, is 
referred to as an “interval.” Herein, a “wave,” such as a P 
wave, will refer to the curvilinear graph (only) portion of a 
time interval, not including the associated time segment. 
Examples of a “wave”, of a “segment,” and of an “interval” 
are illustrated in FIGS. 3A-3B, 3C and 3D-3G, respectively. 

Standard electrocardiography involves multiple record- 
ings of a PQRST complex, referred to as “leads,” which are 
obtained from a plurality o felectrodes or electrode pairs, 
placed at spaced apart locations on a patient’s body. Unipolar 
and bipolar leads are frequently used in standard electrocar- 
diography, for the following purposes: (i) standard bipolar 
limb leads (I, II, III); (ii) augmented unipolar limb leads 
(aVR, aVL, aVF); and (iii) unipolar chest leads (VI, V2, V3, 
V4, V5 V6). As illustrated in FIG. 4, the corresponding elec- 
trode polarities and locations are set forth in Table I. By 
convention, the right leg is treated as “ground.” 


TABLE I 


Lead Polarities and Locations. 

Lead 

Negative Electrode Voltage 

Positive Electrode Voltage 

I 

right arm 

left arm 

II 

right arm 

left leg 

III 

left arm 

left leg 

aVR 

(left arm + left leg)/2 

right arm 

aVL 

(right arm + left leg)/2 

left arm 

aVF 

(right arm + left arm)/2 

left leg 


The bipolar lead voltages are recorded with reference to a 
“ground” electrode located on the right leg. The standard limb 
leads are configured as an equilateral triangle, referred to as 
Einthoven’s triangle, where the following constraint is 
imposed: sum of the voltages impressed for the lead pairs I 
and III is equal to the voltage impressed for the lead pair II. As 
an example, if the QRS impressed voltages for lead pairs I and 
III are 0.8 mV and -0.3 mV, the QRS impressed voltage for 
lead pair II is the algebraic sum, 0.5 mV. 

The augmented unipolar lead voltages are recorded 
between a positive electrode, located on one limb (right arm, 
left arm or left leg), and two negative electrodes, connected 
together and located on the other two limbs (left arm/left leg, 
right arm/left leg and right arm/left arm), respectively. 

The chest lead voltages are recorded between a positive 
electrode located on the patient’s chest and a negative elec- 
trode represented as a sum of voltages for the three standard 
limb electrodes (right arm, left arm, left leg). A sum of the 
three standard limb electrode voltages provides a reference 
value, sometimes referred to as an “indifferent voltage.” The 
locations for the six chest leads are well established in the 
medical field. 

FIG. 4 illustrates placement of some of the electrodes used 
to measure signals and time intervals that are part of an ECG, 
indicating placement of standard limb electrodes on the right 
arm RA, on the left arm LA and on the left leg LL. FIG. 5 
illustrates use of an Einthoven triangle to estimate an atrial 
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depolarization angle 0(ad;depol) associated with a P wave. 
One begins with an equilateral triangle EQTR, with right arm 
(RA), left arm (LA) and left leg (LL) voltages assigned to the 
three vertices as shown. An augmented voltage aVF is mea- 
sured, directed perpendicular to a line segment connecting the 
vertices RA and LA; the vector length of aVF is a deviation of 
the RA measured voltage from the (expected) median value, 
( V (RA)+ V (L A)/2 , described as aVF=V(LL)-{V(RA)+V 
(LA)}/2 

A second augmented voltage aVL is measured perpendicu- 
lar to a line segment connecting the vertices RA and LL, with 
a length represented by a deviation of the measured LA volt- 
age from the expected median value, (V(RA)+V(LL)/2), 
described as aVL=V(LA)-{V(RA)+V(LL)}/2. The atrial 
depolarization vector V(ad;depol) is the vector sum of the 
aVF vector and the aVL vector and is shown in FIG. 5 relative 
to a centroid CT of the triangle EQTR. 

A P-Q time interval, normally of length At(p-q)« 120-200 
msec, represents conduction time from initiation of atrial 
depolarization until initiation of ventricular depolarization, 
which is conventionally measured from the start of the P wave 
to the start of the Q swave. 

Where the ICS is diseased or is affected by presence of 
Digitalis, the P-Q time interval may lengthen as the pulse rate 
decreases; a prolonged P-Q interval, substantially beyond 
200 msec in length, is often evidence of atrio-ventricular 
block. An abnormally short P-Q interval, substantially below 
120 msec in length, is often associated with hypertension 
and/or with paroxysms of tachycardia. The P-Q interval can 
also be shortened where the impulse originates in the AV 
node, or other atrial locations, rather than in the SA node. 

The QRS time interval, normally of temporal length At(q- 
t)«50-100 msec, represents conduction time from initiation 
of ventricular depolarization until initiation of ventricular 
repolarization, and includes spread of the electrical impulse 
through the ventricular muscle. The P signal is normally 
gently rounded, not pointed or notched, and has a temporal 
length «50-l 10 msec. One or more of the P, Q, R, S and/or T 
peaks may be inverted in FIG. 2, depending upon electrode 
placement. A QRS interval greater than about 120 msec in 
temporal length often indicates ventricular arrhythmia or a 
block of one of the bundles. 

Normally, an S-T segment amplitude is approximately 
equal to a T-P segment amplitude. The amplitude of the S-T 
segment, relative to the baseline (e.g., elevated or depressed), 
and the shape of the T signal are often of interest. The T signal 
is normally rounded and slightly asymmetrical. Presence of a 
sharply pointed or grossly notched T signal may indicate 
presence of myocardial infarction (pointed segment) or of 
pericarditis (notched segment). 

In some subjects, a beat (a single PQRST complex) is 
sometimes missed, as illustrated in FIG. 6. This arises from 
the particular physiology of that subject and has not (yet) been 
shown to arise unambiguously from the presence of high 
stress in that subject. 

In normal ECG practice, ten or more electrodes including 
a ground electrode, are attached to the subject at selected, 
spaced apart locations. A chart of each PQRST complex is 
printed separately on a 1 mmxl mm grid, with 1 mm (hori- 
zontal) representing 40 msec (0.04 sec time increment) and 1 
mm (vertical) representing 0.1 millivolts (mV amplitude). An 
upper limit on amplitude is usually 20-30 mm (2-3 mV). The 
chart is normally recorded at a velocity of 25 mm/sec or, 
alternatively, at 50 mm/sec. Measurement of an amplitude of 
5 mm (0.5 mV) or less for all components in a PQRST 
complex is often seen in coronary disease, cardiac failure, 
emphysema and/or obesity. A T signal with unusually large 
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peak amplitude (above 1 mV) may indicate presence of 
ischemia without infarction, or hyperkalemia, or psychosis. 

ECG analysis is generally limited to medical diagnostics 
and to comparison of shifts with the passage of time of ECG 
5 parameters of interest. No substantial work has been done 
applying the ECG results to other areas of interest, such as 
authentication of an asserted identity of a candidate person, 
through analysis of selected ECG results to provide one or 
more physiologically based biometric indicia associated with 
to the candidate person. Nor has any substantial use been made 
of evidence of a malady such as myocardial infarction or 
pericarditis as a characteristic for verifying the identity of a 
candidate person. 

What is needed is a method and associated system for 
1 5 measuring one or more (preferably several) statistical param- 
eters associated with PQRST complexes for a candidate per- 
son and authenticating, or declining to authenticate, the per- 
son’s asserted identity with a reference person, using a 
comparison of the measured statistical parameter values (bio- 
20 metric indicia) with corresponding reference parameter val- 
ues. Optionally, the comparisons should be cumulative so that 
the biometric indicia test can be made progressively more 
demanding, to minimize likelihood of commission of a type I 
error (positive result is false) and/or to balance the likelihoods 
25 of commission of a type I error and commission of a type II 
error (negative result is false) in these comparisons. These 
comparisons should also allow for possible changes with 
passage of time of PQRST complex characteristics for a 
candidate person. Preferably, evidence of presence of a 
30 malady in a reference person should be available for biomet- 
ric use in comparison of a candidate person with this refer- 
ence person. 

SUMMARY OF THE INVENTION 

35 

These needs are met by the invention, which provides a 
method and associated system for authenticating, or declin- 
ing to authenticate, an identity asserted by a candidate person. 
The PQRST complexes for two different persons are likely 
40 sufficiently different that a comparison will allow discrimi- 
nation between these persons for at least one parameter. Fur- 
ther, if a reference person is known to have a cardiac -related 
anomaly, absence of this anomaly in the measured PQRST 
complex of the candidate person strongly indicates that the 
45 candidate person is not the reference person, and conversely; 
presence of this anomaly in the reference person and in the 
candidate person supports, but does not necessarily require, a 
conclusion that the candidate person is the reference person. 

In a first contribution, a sequence of measured values is 
50 provided for at least one peak signal amplitude value, drawn 
from a set of measured peak signal amplitude values, a=ap, 
aq, ar, as and at (positive, negative or zero), for a sequence of 
PQRST heart electrical signal complexes for a candidate per- 
son. A selected combination of one or more statistical param- 
55 eters, drawn from a set of measured parameters, sp 1 (a;meas) 
={normalized mean m„(a;meas)=m(a;meas)/(0.1 mV), 
normalized standard deviation sd„(a;meas)=sd(a;meas)/lm 
(a;meas)l, skewness sk(a;meas), kurtosis ku(a;meas)}, is 
computed, for at least one of the measured peak signal ampli- 
60 tudes , a=ap , aq, ar, as and/ or at . A corresponding combination 
of one or more reference statistical parameters, drawn from a 
set of known reference parameters, spl(a;ref)={normalized 
mean m w (a;ref), normalized standard deviation sd„(a;ref), 
skewness sk(a;ref) and kurtosis ku(a;ref)}, is provided, for 
65 the corresponding reference peak signal amplitude(s), a=ap, 
aq, ar, as and at, for PQRST complexes associated with a 
known reference person. 
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K1 combinations, numbered kl=l, . . . , K1 (Kl^l) of 
selected non-negative weight values, wl, w2, w3 and w4 
(= 1 -wl-w2-w3) are associated with magnitudes of the 
respective differences, m w (a;meas)-m„(a;ref), sd„(a;meas)- 
sd M (a;ref), sk(a.meas)-sk(a;ref) and ku(a;meas_-ku(a;ref), 
where wl , w2, w3 and w4 may depend upon the index kl . A 
first difference Al(kl), dependent upon these differences 
with their associated weights, is computed and compared 
with a first selected range R1 (kl ) of values. When A 1 (kl ) lies 
in the first range Rl(kl) of values, this range condition is 
interpreted as indicating that the candidate person is likely to 
be the reference person; when A1 (kl) does not lie in the first 
range R1 (kl), this condition may be interpreted as indicating 
that the candidate person is not likely to be the reference 
person. 

In a second contribution, the peak signal amplitudes, a=ap, 
aq, ar, as and at, are replaced by time interval lengths, A =At 
(p-q), At(q-r), At(r-s), At(s-t) and At(t-p). A selected combi- 
nation of one or more statistical parameters, drawn from a set 
of measured parameters, sp2(At;meas)={normalized mean 
m w (At;meas)=m(a;meas)/(0.04 sec), normalized standard 
deviation sd w (At;meas)=sd(At;meas)/m(At;meas), skewness 
sk(At;meas), kurtosis ku(At;meas)}, is computed, for at least 
one of the measured time intervals A=At(p-q), At(q-r), At(r-s), 
At(s-t) and At(t-p). The quantity At(t-p) is a time increment 
associated with a refractory period between an end of a ven- 
tricular repolarization signal and a beginning of an atrial 
depolarization signal for the following heart cycle. A corre- 
sponding combination of one or more reference statistical 
parameters, drawn from a set of known reference parameters, 
sp2(At;ref)={normalized mean m M (At;ref), normalized stan- 
dard deviation sd„(At;ref), skewness sk(At;ref) and kurtosis 
ku(At;ref)}, is provided, for the corresponding reference time 
intervals At=At(p-q), At(q-r), At(r-s), At(s-t) and At(t-p), for 
PQRST complexes associated with a known reference per- 
son. 

K2 combinations, numbered k2=l, . . . , K2 (K2^1) of 
selected non-negative weight values, wl’, w2 f , w3’ and 
w4' ( 1 -wl , -w'2-w3') are associated with magnitudes 
of the respective differences, m w (At;meas)-m w (At;ref), 
sd„(At;meas)-sd w (At;ref), sk(At.meas)-sk(At;ref) and ku(At; 
meas)_-ku(At;ref), where wl', w2', w3' and w4' may depend 
upon the index k2. A second difference A2(k2), dependent 
upon these differences with their associated weights, is com- 
puted and compared with a second selected range R2(k2) of 
values. When A2(k2) lies in the second range R2(k2) of 
values, this range condition is interpreted as indicating that 
the candidate person is likely to be the reference person; when 
A2(k2) does not lie in the second range R2(k2), this condition 
may be interpreted as indicating that the candidate person is 
not likely to be the reference person. 

In a third contribution, the measured peak signal amplitude 
values, a=ap, aq, ar, as and at are replaced by depolarization 
and repolarization angles, 0=0 (mode ;ad;depol), 0(mode;vd; 
depol) and 0(mode;vr;repol) (mode=meas or ref), associated 
with the angle 0 that the respective atrial depolarization vec- 
tor V(mode; ad; depol) the ventricular depolarization vector 
V(mode;vd;depol) and the ventricular repolarization vector 
V(mode;vr;repol) makes with a reference baseplane, such as 
a locally horizontal plane, for a sequence of PQRST heart 
electrical signal complexes for a candidate person. A selected 
combination of one or more statistical parameters, drawn 
from a set of measured parameters, sp3(0;meas)={normal- 
ized mean m w (0 ;meas)=m(0 ;meas)/ ( 1 °), normalized standard 
deviation sd w (0;meas)=sd(0;meas)/lm(0;meas)l, skewness 
sk(0;meas), kurtosis ku(0;meas)}, is computed, for at least 
one of the measured angles, 0=0(meas;ad;depol), 0(meas;vd; 
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depol) and 0(meas;vr;repol). A corresponding combination 
of one or more reference statistical parameters, drawn from a 
set of known reference parameters, sp3(0;ref)={normalized 
mean m„(0;ref), normalized standard deviation sd„(0;ref), 
5 skewness sk(0;ref) and kurtosis ku(0;ref)}, is provided, for 
the corresponding reference angles 0=0(ref;ad;depol), 0(ref; 
vd;depol) and 0(ref;vr;repol), for PQRST complexes associ- 
ated with a known reference person. 

K3 combinations, numbered k3=l, . . . , K3 (K3^1) of 
to selected non-negative weight values, wl", w", w3 M and w M 4 
(=1-w m -w2"-w3") are associated with magnitudes of the 
respective differences, m„(0;meas)-m„(0;ref), sd„(0;meas)- 
sd„(0;ref), sk(0.meas)-sk(0;ref) and ku(0;meas_-ku(0;ref), 
where w 1 ”, w2 M , w3 " and w4 M may depend upon the index k3 . 
15 A third difference A3(k3), dependent upon these differences 
with their associated weights, is computed and compared 
with a third selected range R3(k3) of values. When A3(k3) 
lies in the third range R3(k3) of values, this range condition is 
interpreted as indicating that the candidate person is likely to 
20 be the reference person; when A3(k3) does not lie in the first 
range R3(k3), this condition may be interpreted as indicating 
that the candidate person is not likely to be the reference 
person. 

Vectors associated with heart depolarization and repolar- 
25 ization waves are largely predictable, with values that are 
usually close to their mean values. 

In a fourth contribution, the measured peak signal ampli- 
tude values, a=ap, aq, ar, as and at are replaced by a length 
L(mode;ad;depol) of the atrial depolarization vector V(mode; 
30 ad;depol), a length L(mode;vd;depol) of the ventricular depo- 
larization vector V(mode;vd;depol) and a length L(mode;vr; 
repol) of the ventricular repolarization vector V(mode;vr; 
repol) for a sequence of PQRST heart electrical signal 
complexes for a candidate person. A selected combination of 
35 one or more statistical parameters, drawn from a set of mea- 
sured parameters, sp4(L;meas)={normalized mean m w (L; 
meas)=m(L;meas)/(L0), normalized standard deviation sd„ 
(L;meas)=sd(L;meas)/lm(L;meas)l, skewness sk(L;meas), 
kurtosis ku(L;meas)}, is computed, for at least one of the 
40 measured lengths, L=L(meas;ad;depol), L(meas;vd; depol) 
and L(meas;vr;repol). L0 is a reference length of a selected 
one of the vectors V (ref;ad;depol), V (ref ;vd; depol) and V(ref; 
vr;repol), in appropriate units. A corresponding combination 
of one or more reference statistical parameters, drawn from a 
45 set of known reference parameters, sp4(L;ref)={normalized 
mean m„(L;ref), normalized standard deviation sd„(L;ref), 
skewness sk(L;ref) and kurtosis ku(L;ref)}, is provided, for 
the corresponding reference vector lengths L=L(ref;ad;de- 
pol), L(ref;vd;depol) and L(ref;vr;repol), for PQRST com- 
50 plexes associated with a known reference person. 

K4 combinations, numbered k4=l, . . . , K4 (K4i^l) of 
selected non-negative weight values, w ,M l, w2 m , w3 ,M and 
w4 ,M (=l-wl”’-w2 ,M -w3 ,M ) are associated with magnitudes 
of the respective differences, m„(L;meas)-m w (L;ref), sd w (L; 
55 meas)-sd„(L;ref), sk(L.meas)-sk(L;ref) and ku(L;meas_-ku 
(L;ref), where wl ,M , w2 ,M , w3 ,M and w4 ,M may depend upon the 
index k4. A fourth difference A4(k4), dependent upon these 
differences with their associated weights, is computed and 
compared with a fourth selected range R4(k4) of values. 
60 When A4(k4) lies in the fourth range R4(k4) of values, this 
range condition is interpreted as indicating that the candidate 
person is likely to be the reference person; when A1 (k4) does 
not lie in the fourth range R4(k4), this condition may be 
interpreted as indicating that the candidate person is not likely 
65 to be the reference person. 

The difference tests for A1 (kl), A2(k2), A3(k3) and A4(k4) 
may be performed individually and/or may be combined in 
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pairs and/or combined as a triple test to determine if the 
candidate person is likely to be, or is not likely to be, the same 
as the reference person. 

The preceding analysis extends to cross-correlation vari- 
ables <al-a2>, with al, a2={ap, aq, ar, as or at} (al*a2), to 5 
<AtTA2>, with Atl, At2={At(p-q), At(q-r), At(r-s), At(s-t) or 
At(t-p) } (Atl*At2), to <0 1 102> with 01, 03={0(mode;ad; 
depol), 0(mode;vd;depol) or 0(mode;vr;repol)} (01*03), to 
<LTL4>, with LI, L4=L(mode;ad;depol), L(mode;vd;de- 
pol) and L(mode;vr;repol) (L1*L4), and to cross -deviation 10 
variables <al*At2>, <al*03>, <At2*03>, <al-L4>, <At2*03>, 
<At2-L4> and <03*L4>. 

Where an ECG of the reference person indicates presence 
of a non-trivial cardiovascular malady, such as a previous 
myocardial infarction or pericarditis, presence or absence of 15 
this malady in the ECG results of the candidate person can 
provide support for, or refutation of, congruence of the can- 
didate person and the reference person. 

BRIEF DESCRIPTION OF THE DRAWINGS 20 

FIG. 1 schematically illustrates a human heart. 

FIG. 2 schematically illustrates a representative PQRST 
electrical signal complex generated in association with a 
human heart. 25 

FIGS. 3A-3G illustrate “waves”, a “segment” and “inter- 
vals.” 

FIG. 4 schematically illustrates electrode placement for an 
ECG. 

FIG. 5 illustrates use of an Einthoven triangle. 30 

FIG. 6 illustrates a sequence of PQRST complexes, with 
one missed PQRST complex. 

FIGS. 7 A and 7B illustrate positive and negative skewness. 

FIGS. 8, 9, 10 and 11 illustrate procedures to practice the 
invention. 35 

FIG. 12 illustrates the depolarization and repolarization 
vectors, V(mode;ad; depol), V(mode;vd; depol) and V(mode; 
vr;repol) associated with a PQRST complex. 

FIGS. 13 A, 13B and 13C graphically illustrate a method 
for displaying and comparing the measured and reference 40 
parameter values. 

DESCRIPTION OF BEST MODE OF THE 
INVENTION 

45 

In one graphical presentation, a PQRST electrical signal 
curve produced in a cycle of the heart, illustrated schemati- 
cally in FIG. 2, has three maxima and two minima, spaced 
apart in time from each other and has at least 1 6 independent 
associated values: five extremum or peak amplitude values 50 
(ap, aq, ar, as, at), plus five signal-to -signal temporal separa- 
tion values (e.g., At(p-q), At(q-r), At(r-s), At(s-t), At(t-p)), 
measured between time points for two consecutive signals (P, 

Q, R, S, T) of a single PQRST complex or between two 
consecutive complexes, three depolarization and repolariza- 55 
tion angles, 0 (ad; depol), 0(vd;depol) and 0(vr;repol), for 
atrial depolarization (ad) vectors, ventricular depolarization 
(vd) vectors and ventricular repolarization (vr) vectors mea- 
sured relative to a selected baseplane, such as a locally hori- 
zontal plane or locally vertical plane, plus three vector lengths 60 
(magnitudes) L(ad;depol), L(vd; depol) and L(vr;;repol), for 
the ad, vd and vr vectors. Statistical parameters (normalized 
mean, normalized standard deviation, skewness, kurtosis) of 
any subset, or of the full set, of these 1 6 values can be used to 
partly or wholly characterize a person. Because an ECG 65 
works with as many as 12 signal channels, processed from 
signals collected at 10 or more spaced apart electrodes 
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located on a person’s body, as many as 12*16=192 param- 
eters are available for statistical comparison of a candidate 
person and a reference person. 

Consider, for example, the variable aq, which is the curve 
extremum value associated with the q-component of the 
PQRST complex for a particular channel no. b in FIG. 2 
(b=l, . . . , 12), and consider a sequence of measurements of 
peak values aq(j;b) (j=l, 2, . . . , J; J^4) of aq for J heart cycles 
(J^4), not necessarily consecutive. The normalized mean, 
normalized standard deviation, skewness and kurtosis for the 
variable a=ap, aq, ar, as and/or at are defined, respectively, as 

J (l) 

m(a;meas) = I,aq;(j;meas) / J , 

7 = 1 


w„(a;meas)=m(a;meas)/(0.1 mV), (2) 


J (3) 

sd{a\measf‘ — H(a(j;meas ) — m(a\meas)) 2 j {J — 1), 

7=1 


sd n ( a ;meas)=.s<7 ( a ;meas )/m(a ;meas ), (4) 


J (5) 

sk{a\meas) = Y.{a{j\meas) — m{a\meas)y > j (J — \)sd{a\measf , 

7=1 


J (6) 

ku{a;meas ) = Z(a(j;meas) - m{a\meas)f j (J - 1 )sd(a;meas) 4 - 3, 

7=1 

The skewnness coefficient sk(a;meas) is a measure of non- 
symmetry of the variable a about the mean value m(a;meas). 
The skewness for a normally distributed variable is 0. A 
positive value of skewness indicates that the positive branch 
(a>m(a;meas)) of the statistical distribution has a longer tail 
than the negative branch (a<m(a;meas)), and inversely for a 
negative value of skewness, as illustrated in FIGS. 7 A and 7B. 

The kurtosis coefficient ku(a;meas) is a measure of “sharp- 
ness of peak” of the variable a, when compared to a normally 
distributed variable, which has a kurtosis of 0. A positive (or 
negative) value of kurtosis indicates that the variable a is more 
sharply peaked (or less sharply peaked) than a corresponding 
normally distributed variable. 

In a first contribution, a sequence of measured values is 
provided for at least one peak signal amplitude value, drawn 
from a set of measured peak signal amplitude values, a=ap, 
aq, ar, as and at (positive, negative or zero), for a sequence of 
PQRST heart electrical signal complexes for a candidate per- 
son. A selected combination of one or more statistical 
parameters, drawn from a set of measured parameters, spl (a; 
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meas)={nonnalized mean m w (a;meas)=m(a;meas)/(0. 1 mV), 
normalized standard deviation sd w (a;meas)=sd(a;meas)/lm 
(a;meas)l, skewness sk(a;meas), kurtosis ku(a;meas)}, is 
computed, for at least one of the measured peak signal ampli- 
tudes, a=ap, aq, ar, as and/or at. A corresponding combination 
of one or more reference statistical parameters, drawn from a 
set of known reference parameters, spl(a;ref)={normalized 
mean m„(a;ref), normalized standard deviation sd„(a;ref), 
skewness sk(a;ref) and kurtosis ku(a;ref)}, is provided, for 
the corresponding reference peak signal amplitude(s), a=ap, 
aq, ar, as and at, for PQRST complexes associated with a 
known reference person. 

K1 combinations, numbered kl=l, . . . , K1 (Kli^l) of 
selected non-negative weight values, wl, w2, w3 and w4 
(=l-wl-w2-w3) are associated with magnitudes of the 
respective differences, m M (a;meas)-m w (a;ref), sd w (a;meas)- 
sd M (a;ref), sk(a.meas)-sk(a;ref) and ku(a;meas_-ku(a;ref), 
where wl , w2, w3 and w4 may depend upon the index kl . A 
first difference Al(kl), dependent upon these differences 
with their associated weights, is computed and compared 
with a first selected range R1 (kl ) of values. When A 1 (kl ) lies 
in the first range Rl(kl) of values, this range condition is 
interpreted as indicating that the candidate person is likely to 
be the reference person; when A1 (kl) does not lie in the first 
range R1 (kl), this condition may be interpreted as indicating 
that the candidate person is not likely to be the reference 
person. 

The system also measures the variables At=At(p-q;meas), 
At(q-r;meas), At(r-s;meas), At(s-t;meas) and/or At(t-p;meas), 
The quantity At(t-p) is a time increment associated with a 
refractory period between an end of a ventricular repolariza- 
tion signal and a beginning of an atrial depolarization 
signal for the following heart cycle. The system computes one 
or more of the corresponding statistical parameters m w (At; 
meas)=m(At;meas)/ (0.04 sec), sd w (At;meas)=sd(At;meas)/m 
(At;meas), sk(At;meas) andku(At;meas), drawn from a set of 
statistical parameters sp2(At;meas)={m„(At;meas), sd„(At; 
meas), sk(At;meas), ku(At;meas)}, and provides the corre- 
sponding statistical parameters m„(At;ref)=m(At;ref)/(0.04 
sec), sd w (At;ref)=sd(At;ref)/m(At;ref), sk(At;ref) and 
ku(At;ref), drawn from a set of statistical parameters sp2(At; 
ref)={m w (At;ref), sd w (At;ref), sk(At;ref), ku(At;ref)}. 

K2 combinations, numbered k2=l, . . . , K2 (Kli^l) of 
selected non-negative weight values, wl', w2’, w3' and 
w4' (=l-wT-w'2-w3') are associated magnitudes of 
with the respective differences, m„(At;meas)-m w (At;ref), sd„ 
(At;meas)-sd w (At;ref), lsk(At.meas)-sk(At;ref) and ku(At; 
meas)_-ku(At;ref), where wl', w2', w3' and w4' may depend 
upon the index k2. 

A second difference A2(k2), dependent upon these differ- 
ences with their associated weights, is computed and com- 
pared with a second selected range R2(k2) of values. When 
A2(k2) lies in the second range R2(k2) of values, this range 
condition is interpreted as indicating that the candidate per- 
son is likely to be the reference person; when A2(k2) does not 
lie in the second range R2(k2), this condition may be inter- 
preted as indicating that the candidate person is not likely to 
be the reference person. 

The system also measures the angle variables 0(mode;ad; 
depol), 0(mode;vd;depol) and 0(mode;vr;repol), associated 
with the angle 0 that the respective atrial depolarization vec- 
tor V(mode;ad; depol), the ventricular depolarization vector 
V(mode;vd;depol) and the ventricular repolarization vector 
V(mode;vr;repol) makes, respectively, with a reference base- 
plane, such as a locally horizontal plane. The system com- 
putes one or more statistical parameters, drawn from a set of 
measured parameters, sp3 (0 ;meas)={normalized mean m„(0 ; 
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meas)=m(0;meas)/(l °), normalized standard deviation sd„(0; 
meas)=sd(0;meas)/lm(0;meas)0, skewness sk(0;meas), kur- 
tosis ku(0;meas)} and provides corresponding statistical 
parameters, drawn from a set of reference parameters sp3(0; 
5 ref)={normalized mean m w (0 ;ref)=m(0 ;ref)/ ( 1 °), normalized 
standard deviation sd„(0;ref)=sd(0;ref)/lm(0;ref)l, skewness 
sk(0;ref), kurtosis ku(0;ref)}. 

K3 combinations, numbered k3=l, . . . , K3 (K3^1) of 
selected non-negative weight values, wl", w", w3 M and w M 4 
10 (=1-w m -w2"-w3") are associated with magnitudes of the 
respective differences, m w (0;meas)-m w (0;ref), lsd w (0;meas)- 
sd„(0;ref), sk(0.meas)-sk(0;ref)l andku(0;meas_-ku(0;ref), 
where w 1", w2 M , w3 " and w4 M may depend upon the index k3 . 
15 A third difference A3(k3), dependent upon these differ- 
ences with their associated weights, is computed and com- 
pared with a second selected range R3(k3) of values. When 
A3(k3) lies in the third range R3(k3) of values, this range 
condition is interpreted as indicating that the candidate per- 
20 son is likely to be the reference person; when A3 (k3) does not 
lie in the third range R3(k3), this condition may be interpreted 
as indicating that the candidate person is not likely to be the 
reference person. 

The system also measures the vector length variables 
25 L(mode;ad;depol), L(mode;vd;depol) and L(mode;vr;repol), 
associated with the atrial depolarization vector V(mode;ad; 
depol), the ventricular depolarization vector V(mode;vd;de- 
pol) and the ventricular repolarization vector V(mode;vr;re- 
pol) respectively. A selected combination of one or more 
30 statistical parameters, drawn from a set of measured param- 
eters, sp4(L;meas)={normalized mean m„(L;meas)=m(L; 
meas)/(L0), normalized standard deviation sd„(L;meas)=sd 
(L;meas)/lm(L;meas)l, skewness sk(L;meas), kurtosis ku(L; 
meas)}, is computed, for at least one of the measured lengths, 
35 L=L(meas;ad;depol), L(meas;vd;depol) and L(meas;vr;re- 
pol). L0 is a reference length of a selected one of the vectors 
V(ref;ad;depol), V(ref;vd;depol) and V(ref;vr;repol), in 
appropriate units. A corresponding combination of one or 
more reference statistical parameters, drawn from a set of 
40 known reference parameters, sp4(L;ref)={normalized mean 
m„(L;ref), normalized standard deviation sd„(L;ref), skew- 
ness sk(L;ref) and kurtosis ku(L;ref)}, is provided, for the 
corresponding reference lengths L=L(ref;ad;depol), L(ref; 
vd; depol) and L(ref;vr;repol), for PQRST complexes associ- 
45 ated with a known reference person. 

K4 combinations, numbered k4=l, . . . , K4 (K4i^l) of 
selected non-negative weight values, wl ,M , w ,M , w3 ,M and w M, 4 
(=1-w ,m -w2 ,m -w 3 ,m ) are associated with magnitudes of the 
respective differences, m ( (L;meas)-m w (L;ref), sd w (L;meas)- 
50 sd„(L;ref), sk(L.meas)-sk(L;ref) andku(L;meas_-ku(L;ref), 
where wl ,M , w2 ,M , w3 ,M and w4 ,M may depend upon the index 
k4. 

A fourth difference, A4(k4), dependent upon magnitudes 
of one or more of the differences, m w (L;meas)-m„(L;ref), 
55 sd„(L;meas)-sd„(L;ref), sk(L;meas)-sk(L;ref), and/or ku(L; 
meas)-ku(L;ref), is computed and compared with a fourth 
selected range R4(k4) of values. When A4(k4) lies in the 
fourth range R4(k4) of values, this condition is interpreted as 
indicating that the candidate person is likely to be the refer- 
60 ence person; when A4(k4) does not lie in the fourth range 
R4(k4), this condition may be interpreted as indicating that 
the candidate person is not likely to be the reference person. 
The values Al(kl), A2(k2), A3(k3) and/or A4(k4), relative to 
the respective ranges Rl(kl), R2(k2), R3(k3) and/or R4(k4), 
65 can be displayed, graphically or alphanumerically, or another 
indication can be displayed indicating that the candidate per- 
son is likely to be, or is not likely to be, the reference person. 
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The functional forms of the first, second, third and fourth 
differences, Al(kl), A2(k2), A3(k3) and A4(k4), may be cho- 
sen independently of each other, as functions of differences of 
the corresponding statistical measured and statistical refer- 
ence values. One general class of difference functions is the 
Minkowski distance, defined as follows. 

A1 (kl)={w\ lm„(a;meas)-m„(fl;ref) \^+w2 1 sd„(a; 

meas)-sd n (a;ref) l^+w3 U£(fl;meas)-s£(fl;ref)K 1 + 
w4 lku(a;meas)-ku(a;ref) (7) 

A2(k2)={wl'lm n (At;meas)-m n (At;ref)l /2 +w2'lsd„(At; 
meas)-sd„(At;ref) w3 1 U&(A?;meas)-.s , A:(Af;:ref) 

Z 2 +w4 1 \ku(At;mea.s)-ku(At;ref) Z 2 } 1/ - /2 , (8) 

A3(^3)={wr , lm„(0;meas)-m„(0;ref)^ 3 w2 ,, U<i„(0; 

meas)-5^„(0;ref) I /3 +w 3" UA:(0;meas)-5A:(0;ref) l /3 + 
w4"lfe/(0;meas)-^w(0;ref)K 3 } 1/ - /3 , (9) 

A4(k4)={w\ \m n (L; meas) - m n (L ;ref) f x +w2'" \sd n (L ; 
meas)-5i7„(Z;ref) K 4 +w3 ,M I sk(L ; meas)-.s&(Z ;ref) 
l^+w4 ,,, l^w(Z;meas)-^w(Z,;ref)l /2t } 1/ '^, (10) 

where fl, f2, f3 and f4 are independently chosen positive 
numbers. In Eq. (7) one can prove that, if a particular term, 
wm*lf m (a;rneas)-f m (a;ref)l, has wm>0 and the multiplier lf m 
(a;meas)-f m (a;ref)l is larger than any other multiplier in that 
sum, the difference Al(kl) tends to the value ||f w (a;meas)-f m 
(a;ref)l as the power fl increases without limit (fl— »oo). An 
analogous result is obtained for the difference A2(k2), for the 
difference A3(k3), and for the difference A4(k4). A choice 
fl =2 or f2=2 or f3=2 or f4=2 reproduces a weighted Euclid- 
ean distance in Eq. (7), Eq. (8), Eq. (9) or Eq. (10). Other 
functional forms for Al(kl), A2(k2), A3(k3) and A4(k4) can 
also be chosen. A choice fl=l or f!2=l or f3=l or f4=l repro- 
duces a weighted quasilinear distance. Use of a maximum 
likelihood estimator to estimate a log-likelihood statistical 
distance for verification is discussed in an Appendix. 

Where two or more peak signal amplitudes, for example, 
a=ap and a=ar, are separately examined, the differences 
A 1 (kl ) are computed separately for each signal amplitude (ap 
and ar), using a selected functional form for these differences. 
For the five- fold set of peak signal amplitudes (ap and/or aq 
and/or ar and/or as and/or at), the number of different com- 
binations available is 2 5 - 1=31. 

Where two or more time intervals, for example, At=At(q-r) 
and At=At(t-p), are separately examined, the differences 
A2(k2) are computed separately for each signal amplitude 
(At(q-r) and At(t-p)), using a selected functional form for 
these differences. For the five- fold set of time intervals 
(At(p-r) and/or At(q-r) and/or At(r-s) and/or At(s-t) and/or 
At(t-p)), the number of different combinations available is 
also 2 5 -l=31. 

Where two or more depol/repol angles, for example, 0=0 
(mode;ad;depol) and 0=0 (mode ;vr;repol), are separately 
examined, the differences A3(k3) are computed separately for 
each angle, 0(mode;ad;depol) and 0(mode;vr;repol), using a 
selected functional form for these differences. For the three- 
fold set of angles (0(mode;ad;depol) and/or (0(mode;vd;de- 
pol) and/or 0(mode;vr;repol), the number of different com- 
binations available is 2 3 -l=7. In a similar manner, for the 
three-fold set of vector lengths L(mode;ad;depol) and/or 
L(mode;vd;depol) and/or L(mode;vr;repol), the number of 
available combinations is 2 3 -l=7. 

In a first embodiment, where a first difference Al(kl) lies 
within a first selected range Rl(kl) of difference values for 
each of these Kl combinations; or a second difference A2(k2) 
lies within a second selected range R2(k2) of difference val- 
ues for each of these K2 combinations; or a third difference, 
A3(k3) lies within a third selected range R3(k3) of difference 
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values for each of these K3 combinations, or a fourth differ- 
ence, A4(k4) lies within a fourth selected range R4(k4) of 
difference values for each of these K4 combinations (referred 
to as “range conditions”), the candidate person is determined 
5 to be likely to be the reference person. Satisfaction of at least 
one of the first range condition for Al(kl), the second range 
condition for A2(k2), the third range condition for A3(k3), 
and the fourth range condition for A4(k4) is interpreted as 
indicating that the candidate person and the reference person 
10 are the same person. Where Al(kl) does not lie in the first 
range R1 (kl ), A2(k2) does not lie in the second range R2(k2), 
A3 (k3) does not lie in the third range R3 (k3), and A4(k4) does 
not lie in the fourth range R4(k4) this condition may be 
1 5 interpreted as indicating that the candidate person is likely not 
the reference person. 

In a second embodiment, where at least two of the four 
range conditions on Al(kl), A2(k2), A3(k3) and A4(k4) are 
satisfied, these conditions are interpreted as indicating that 
20 the candidate person and the reference person are likely the 
same person. Where no more than one of the four range 
conditions on A1 (kl), A2(k2), A3(k3) and A4(k4) is satisfied, 
this may be interpreted as indicating that the candidate person 
is likely not the reference person. 

25 In a third embodiment, where at least three of the four range 

conditions on Al(kl), A2(k2), A3(k3) and A4(k4) are satis- 
fied, these conditions are interpreted as indicating that the 
candidate person and the reference person are likely the same 
person. Where no more than two of the four range conditions 
30 on Al(kl), A2(k2), A3(k3) and A4(k4) are satisfied, this may 
be interpreted as indicating that the candidate person is likely 
not the reference person. 

In a fourth embodiment, where all four of the range con- 
ditions on Al(kl), A2(k2), A3(k3) and A4(k4) are satisfied, 
35 these conditions are interpreted as indicating that the candi- 
date person and the reference person are likely the same 
person. Where no more than three of the four range conditions 
on Al(kl), A2(k2), A3(k3) and A4(k4) are satisfied, this may 
be interpreted as indicating that the candidate person is likely 
40 not the reference person. 

For a given (known) person, it may be appropriate to pro- 
vide two or more combinations of references parameters, 
spl(a;ref) and/or two or more sets of reference parameters 
sp2(At;ref) and/or two or more sets of reference parameters 
45 sp3(0;ref) and/or two or more sets of reference parameters 
sp4(L;ref), corresponding to different situations in which 
these parameters were determined, for comparison with cor- 
responding measurements for the candidate person, for at 
least two reasons. Where two or more sets of reference person 
50 parameters are available, the measured parameters for the 
candidate person are preferably compared with each of the 
reference person parameter sets to determine if at least one of 
the reference person parameter sets is in substantial agree- 
ment with the measured candidate person set, as reflected in 
55 the difference values Al(kl), A2(k2), A3(k3) and/or A4(k4). 

First, presence of the well known “white coat syndrome” in 
some reference situations but not in others, may cause an 
end-to-end reference time interval length At(p-p) to become 
significantly smaller in one reference situation than in another 
60 situation, for a given reference person. This will cause at least 
one of the reference time interval lengths, At(p-q), At(q-r), 
At(r-s), At(s-t) and At(t-p), to also become smaller, among 
other changes. Thus, it may be appropriate to compare two or 
more measured combinations Cl(kl)(a;meas) and C2(k2) 
65 (At;meas), C3(k3)(0;meas) and C4(k4)(L;meas) against each 
of two or more corresponding reference combinations for 
different situations for a given candidate person. 
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Second, the cardio -physiology of the candidate person 
may change appreciably with time so that comparison of 
presently measured ECG parameters with corresponding ref- 
erence parameters, prepared several months or years earlier, 
may lead to an increased incidence of Type II errors, wherein 5 
the system erroneously declines to identify the candidate 
person with the reference person. 

The locations on the body of the ECG electrodes used to 
measure the ECG plots for the reference parameter values, 
x(ref) and At(ref) and 0(ref), should be noted carefully and 10 
preserved with the reference parameter values themselves. 
The same body locations should be used for measuring the 
PQRST complex parameters for the candidate person, to 
minimize dispersion problems, amplitude distortion prob- 15 
lems and/or phase distortion problems associated with trans- 
mission of electrical signals through portions of, or along 
surfaces of, the human body. 

Cross -correlations of two or more of the peak amplitude 
signals ap, aq, ar, as and at, and of the time interval lengths 20 
At(p-q), At(q-r), At(r-s), At(s-t) and At(t-p), and of the propa- 
gation angles 0(mode;ad;depol), 0(mode;vd;depol) and 
0(mode;vr;repol), and of the vector lengths L(mode;ad;de- 
pol), L(mode;vd;depol) and L(mode;vr;repol), can also pro- 
vide reference parameters, against which the corresponding 25 
measured parameters of the candidate person can be com- 
pared. A first such parameter is the correlation of two distinct 
peak amplitude signals, 

30 

J (11) 


{ al ■ a2 ) = Lal(j)-a2(j)/(J - 1), 
7 = 1 


35 


al, al=ap, aq, ar, as or at (ah*a2). (12) 

A second such parameter is the correlation of two distinct 
time interval lengths, 

J (13) 

<Arl ■ Ar2> = ZArl (j) ■ Af2(y) / ( J - 1), 

7 = 1 


Atl,At2=At(p-q), A t(q-r), A t(r-s), At(s-t) or A t(t-p) 
(Afl*Af2). 


50 

(14) 


A third such parameter is the correlation of two of the three 
propagation angles 


55 


J 


(15) 


< 01 ■ 03 > = 201 (y) ■ 93(j)f(J - 1) 
7 = 1 


60 


01, 03=0(mode;a;depol), 0(mode;v;depol), 0(mode;v; 

repol) (01*03). (16) 

A fourth such parameter is the correlation of two of the 
three vector lengths 


J (17) 

<L1 -13) = 2Al(/)-L3(/)/ (7 - 1) 

7=1 


LI, Z3=Z,(mode;aa?;depol), Z,(mode;vaf;depol), Z,(mn- 

ode;vr;repol) (Z,*Z,3). (18) 

Any of the set of ( 5 2 )=10 correlation parameters <al-a2> 
and/or any of the set of ( 5 2 ) = 10 correlation parameters 
<Atl -At2> and/or any of the set of ( 3 2 )=3 correlation param- 
eters <01*03> and/or any of the set of ( 3 2 )=3 correlation 
parameters <L1 *L4> can be used as a biometric indicium for 
comparison of a measured value (for a candidate person) 
against a corresponding reference value. Again, the differ- 
ence between the measured parameter value and the corre- 
sponding reference value should fall into a selected range, if 
the candidate person is to be identified as the reference per- 
son. Further, one or more of the 26 statistical parameters, 
m„(a), sd„(a), sk(a), ku(a), m„(At), sd„(At), sk(At), ku(At), 
m„(0), s(£(0), sk(0), ku(0), m„(L), sd„(L), sk(L), ku(L), 
<al*a2>, <Atl*At2>, <01*04>, <L1\L2> <al-At2>, <al-03>, 
<al*L4>, <At2*03>, <At2*L4> and <03*L4> can be used for 
comparison purposes. 

Two or more of the peak signal amplitude variables, a=ap, 
aq, ar, as and at, may be correlated within a single PQRST 
complex, because of constraints on signal recovery (positive- 
to-negative or negative-to-positive) or because of correlations 
between signal distortions introduced by signal transmissions 
through one or more body organs or tissues. Two or more of 
the time interval lengths, At=At(p-q), At(q-r), At(r-s), At(s-t) 
and At(t-p), may be correlated (within a single PQRST com- 
plex or between two consecutive PQRST complexes), 
because of an overall time constraint that is present. Two or 
more of the propagation angles, 01=0 (mode ;ad;depol), 02=0 
(mode;vd;depol), and 03=0 (mode ;vr;repol), may be corre- 
lated because of geometrical or material constraints, associ- 
ated with signal propagation adjacent to the heart. Two or 
more of the propagation vector lengths, Ll=L(mode;ad;de- 
pol), L2=L(mode;vd;depol), and L3=L(mode;vr;repol), may 
be correlated because of geometrical or material constraints, 
associated with signal propagation adjacent to the heart. 

The quantity <xl*At2>, with (xl,At2)={(ap,At(p-q)), (aq, 
At(p-q)), (aq,At(q-r)), (ar,At(q-r)), (ar,At(r-s), (as,At(r-s)), (as, 
At(s-t)), (at,At(s-t)), (at,At(t-p)) or (ap,At(t-p))}, may have 
non-zero cross-deviation, because the magnitude of a peak 
(e.g., aq) may be influenced by at least one of the adjacent 
time interval lengths (e.g., At(p-q) and At(q-r)). Similar non- 
zero cross-deviations, <al*03>, <al*L4>, <At2*03>, 

<At2*L4> and <03-L4>, may occur because of interactions of 
the various ECG variables. 

One can also use cross -deviations of peak amplitudes, time 
interval lengths and the propagation angles for comparison 
purposes, such as the 24 parameters 

<al ■At2>='Z J al (j)-At2(j)/(J- 1 ), (19) 

(al,At2)={ap-At(p-q), aq-At(p-q), xq-At(q-r'), ar-At(q- 
r), ar-At(r-s'), as-At(r-s'), as-At(s-t), at-At(s-f), 
at-At(t-p) or ap-At(t-p)} (20) 

<a 1 -3 >=2yfl 1 (/) -03 (f)/ (J- 1 ) , (21) 


(al,03)={fl/>-0(a<a?;depol), a^-0(fl(7;depol), aq-Q(vd; 
depol), flr-0(vi7;depol), A5-0(vi7;depol), as-Q(vr; 
repol), fl?-0(vr;repol)}, 


( 22 ) 
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<a 1 -14 >=2/i 1 (j)IA(j)!(J- 1 ), (23) 

(al^A)={ap-L{ad;A&po\), aq-L(ad; depol), aq-L(vd; 
depol), ar-L(vd; depol), as-L(vd; depol), as-L(vr; 
repol), at-L(vr; repol)}, (24) 

5 

<At2-03>=ZjAt2(j)‘03 (j)/(J- 1), (25) 

(A?2,03)={A?(p-#)-0(flrf;depol), A?(#-r)-0(v^;depol), 

A?(r-5)-0(vrf;depol), Ar(5-?)-0(v^;depol), Af(.s-r)-0 
(vr;repol), A?(*-p)-0(vr;repol), A?(?-p)-0(ai/;de- 

pol)} (26) 10 

<At2-L4->=H j At2(J) ) L(J)/(J-l), (27) 

(Af2,Z4)={A?(p-#)-Z,(flrf;depol), Ar(#-r)-Z(v</;depol), 

At(r-s)-L {yd ; dep o 1) , A?(s-f)’Z(vrf;depol), ^ 

A?(5-?)-Z(vr;repol), A?(?-p)-Z(vr;repol), At(t-p)-L 

(a<7;depol)}, (28) 

<03 ■Z4>=2 / -03 (f)L4(j)l(J- 1 ) (29) 

(03,Z4)={0(flaf;depol)-Z(a^;depol), 0(v^;depol)-Z(a^; 20 

depol), Q(ad; dep ol) -L(vd; dep ol) , 0(v<7;depol)-Z 
(vi7; depol), 0(v^;depol)-Z(vr;repol), 0(vr;repol)-Z 
(vi7; depol), 0(vr;repol)-Z(va?;repol)}. (30) 


One can also use certain related functions of the correla- 
tions for comparison purposes. Consider, for example, a 25 
cross -correlation function 


J (31) 

12(A) = Z(al(/) - X-a2(j)) 2 !(J - 1) 

7 = 1 

= <al 2 > - 2A<al -x2) + X 2 (al 2 ) 


pi 2(A=1 )=<a l 2 >-2 <a 1 -a 2>+<a2 2 >. (32) 

The function p,12(X) is minimized by the choice 

A=A(mm)=<al-a2>/<a2 2 >, (33) 

pi 2(A=A(min))=<a l 2 >-{ <a 1 ■a2>} 2 /<a2 2 >. (34) 

Another quantity of interest is a cross- standard deviation, 
defined as 


J (35) 

var(al , al) = I,(al(j) - m(al )) ■ ( a2(j ) - m(a2))/(J - 1). 


= <al ■ al) - m(al ) ■ m(a2). 

Similar results are obtained where the variable a is replaced 
by the variable At or by the variable 0. 55 

FIG. 8 is a flow chart illustrating a first embodiment of the 
invention. In step 81 , at least one of the difference values, 
Al(kl), A2(k2), A3(k3) and A4(k4) is computed and com- 
pared with the corresponding range, Rl(kl), R2(k2), R3(k3) 
and R4(k4), respectively. The system determines if A 1 (kl ) is 60 
present (computed) and lies in the first range Rl(kl), in step 
82 ; if A2(k2) is present (computed) and lies in the second 
range R2(k2), in step 83 , if A3(k3) is present (computed) and 
lies in the third range R3(k3), in step 84 , and/or if A4(k4) is 
present (computed) and lies in the fourth range R4 (k4), in step 65 
85 . The number M of “yes” answers in steps 82-85 is accu- 
mulated in step 86 . If at least one of the answers to the queries 
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in step 82 , 83, 84 and 85 is “yes” (M^ 1 in step 87 ), the system 
interprets this condition as indicating that the candidate per- 
son is likely the reference person, in step 88. If none of the 
answers to the queries in step 82 , step 83 , step 84 and step 85 
is “yes,” the system optionally interprets this condition as 
indicating that the candidate person is likely not the reference 
person, in step 89 . 

FIG. 9 is a flow chart illustrating a second embodiment of 
the invention. In step 91 , at least two of the difference values, 
Al(kl), A2(k2), A3(k3) and A4(k4) are computed and com- 
pared with the corresponding ranges, R1 (kl ), R2(k2), R3(k3) 
and R4(k4), respectively. The system determines if Al(kl) is 
present (computed) and lies in the first range Rl(kl), in step 
92 ; if A2(k2) is present (computed) and lies in the second 
range R2(k2), in step 93 , if A3(k3) is present (computed) and 
lies in the third range R3(k3), in step 94 , and/or if A4(k4) is 
present (computed) and lies in the fourth range R4(k4), in step 
95 . The number M of “yes” answers in steps 92-95 is accu- 
mulated in step 96 . If at least two of the answers to the queries 
in step 92 , 93, 94 and 85 is “yes” (Mi^2 in step 97 ), the system 
interprets this condition as indicating that the candidate per- 
son is likely the reference person, in step 98 . If zero or one of 
the answers to the queries in step 92 , step 93 , step 94 and step 
95 is “yes,” the system optionally interprets this condition as 
indicating that the candidate person is likely not the reference 
person, in step 99 . 

FIG. 10 is a flow chart illustrating a third embodiment of 
the invention. In step 101 , at least three of the difference 
values, Al (kl), A2(k2), A3(k3) and A4(k4) are computed and 
compared with the corresponding ranges, Rl(kl), R2(k2), 
R3(k3) and R4(k4), respectively. The system determines if 
Al(kl) is present (computed) and lies in the first range 
Rl(kl), in step 102 ; if A2(k2) is present (computed) and lies 
in the second range R2(k2), in step 103 , if A3(k3) is present 
(computed) and lies in the third range R3(k3), in step 104 , 
and/or if A4(k4) is present (computed) and lies in the fourth 
range R4(k4), in step 105 . The number M of “yes” answers in 
steps 102-105 is accumulated in step 106 . If at least three of 
the answers to the queries in step 102 , 103 , 104 and 105 is 
“yes” (Mi^3 in step 107 ), the system interprets this condition 
as indicating that the candidate person is likely the reference 
person, in step 108 . If zero, one or two of the answers to the 
queries in step 102 , step 103 , step 104 and step 105 is “yes,” 
the system optionally interprets this condition as indicating 
that the candidate person is likely not the reference person, in 
step 109 . 

FIG. 11 is a flow chart illustrating a fourth embodiment of 
the invention. In step 111 , at least one of each of the four 
difference values, Al (kl), A2(k2), A3(k3) and A4(k4) is com- 
puted and compared with the corresponding ranges, Rl(kl), 
R2(k2), R3(k3) and R4(k4), respectively. The system deter- 
mines if Al(kl) is present (computed) and lies in the first 
range R1 (kl), in step 112 ; if A2(k2) is present (computed) and 
lies in the second range R2(k2), in step 113 , if A3(k3) is 
present (computed) and lies in the third range R3(k3), in step 
114 , and/or if A4(k4) is present (computed) and lies in the 
fourth range R4(k4), in step 115 . The number M of “yes” 
answers in steps 112-115 is accumulated in step 116 . If each 
of the answers to the queries in step 112 , 113, 114 and 115 is 
“yes” (M^4 in step 117 ), the system interprets this condition 
as indicating that the candidate person is likely the reference 
person, in step 118 . If fewer than four of the answers to the 
queries in step 112 , step 113 , step 114 and step 115 is “yes,” 
the system optionally interprets this condition as indicating 
that the candidate person is likely not the reference person, in 
step 119 . 
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Comparison of each of, or a subgroup of, the 1 5 combina- 
tions, spl(a;meas), sp2(At;meas), sp3(0;meas) and sp4(L; 
meas), with the corresponding combinations, spl(a;ref), sp2 
(At;ref), sp3(0;ref) and sp4(L;ref), allows discrimination 
between a candidate person and each of a large number of 
reference persons. For example, if each of the 1 5 reference 
combinations, spl(a;ref), sp2(At;ref), sp3(0;ref) and sp4(L; 
ref), has six discrete or distinguishable values, the number of 
different combinations is 6 15 , or about 470 billion. Supple- 
menting these reference measures with the reference cross 
correlations and cross deviations, <xl*x2;ref>, <Atl*At2; 
ref>, <0T03;ref>, <LTL2;ref>, <xl-At2;ref>, <xl*03;ref>, 
<xl-L4;ref>, <Atl*03;ref>, <At2*L4> and <03*L4;ref> 
increases the number of combinations by many orders of 
magnitude so that in principle each person in a population 
much larger than 1.2xl0 12 can be distinguished, using these 
measures. 

Where a reference person is known to have a cardiac- 
related anomaly or malady that is not transitory (doesn’t 
disappear or sharply diminish with passage of time), and this 
artifact is manifested by a particular pattern in the PQRST 
complex, absence of this particular pattern in a PQRST com- 
plex presented by a candidate person supports a conclusion 
that the candidate person is not the reference person, irrespec- 
tive of agreement or disagreement of other measured statis- 
tical parameter values and corresponding reference values. A 
similar conclusion may be drawn, but is not required, where 
the candidate person has a non-transitory, cardiac -related 
anomaly or malady and the reference person does not have the 
corresponding anomaly or malady. 

For example: At(p-r;ref)»200 msec is associated with 
atrio -ventricular block; At(p-r;ref)<120 msec is associated 
with hypertension and/or paroxysms of tachycardia; At(q-s; 
ref)>120 msec is associated with ventricular arrhythmia or a 
block of one of the bundles; and presence of a sharply pointed 
or grossly notched t segment is associated with myocardial 
infarction (pointed segment) or with pericarditis (notched 
segment). 

The converse is not necessarily true: presence of a particu- 
lar pattern (anomaly or malady) in a statistical measure, but 
not in the corresponding pattern for a reference person may 
support, but does not necessarily require, a conclusion that the 
candidate person is not the reference person. The reference 
person may have developed this pattern after the most recent 
reference pattern was formed. 

The ten independent parameters, ap, aq, ar, as, at, At(p-q), 
At(q-r), At(r-s), At(s-t) and At(t-p), discussed in the preceding 
belong to a set of about 192 heart signal parameters, not all of 
which may be independent. This includes, for example, a 
post-t plateau region, having an amplitude denoted as a(t; 
post) (FIG. 2). Ideally, the value of a(t;post) is 0 mV. Where 
a(t;post)^-0.5 mV, this condition indicates that the subject 
has recently experienced a myocardial infarction, probably 
within the preceding two weeks. Where a(t;post)^-0.5 mV at 
a given time, the magnitude la(t;post)l will decrease with 
passage of time so that evidence of an earlier myocardial 
infarction will slowly disappear. An amplitude a(t;post) with 
a substantial non-zero magnitude, e.g., Ia(t;post)l^0.5 mV, 
will appear graphically on a only if the frequency range for 
signal processing extends down to about 0.05 Hz, or lower. 
Where the frequency range does not extend below about 1 Hz, 
for example, the amplitude a(t;post) will often have a sub- 
stantially zero magnitude. 

Three other heart signal parameters of interest are the 
direction of a vector V(mode;ad;depol) for sino atrial-to- 
atrial-ventricular depolarization, associated with the p seg- 
ment, the direction of a vector V(mode;vd;depol) for sino 
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atrial -to-atrial -ventricular depolarization, associated with the 
qrs segment, and the direction of a vector V(mode;vr;repol) 
for atrial-ventricular repolarization, associated with the t seg- 
ment, which are expressed as angles 0(mode;ad;depol), 
5 0(mode;vd;depol) and 0(mode;vr;repol), respectively, rela- 
tive to a horizontal baseline for the subject, as illustrated in 
FIG. 2. For a given subject, 0(mode;ad;depol), 0(mode;vd; 
depol) and 0(mode;vr;repol) will have values in ranges of 
0°-90°, and these angles do not vary significantly from one 
to heart cycle to the next, absent a deleterious change in the 
cardio condition of the subject. The angles 0(mode;ad;depol), 
0(mode;vd;depol) and 0(mode;vr;repol) for a given subject 
appear to be independent of each other and can be measured 
to within an inaccuracy of a few degrees. Thus, for example, 
15 an angular range of O°^0^9O° can be decomposed into 6 
central angle values, spaced about 15° apart. At least six 
length values (magnitudes) for each of the corresponding 
vectors, V(ad; depol), V(vd; depol) and V(vr;repol), are esti- 
mated to be distinguishable so that angle and length measure- 
20 ments provide at least 36 distinguishable value pairs. 

This last number (36), multiplied by an estimated number 
of independent multi-values of the set {ap, aq, ar, as, at, 
At(p-q), At(q-r), At(r-s), At(s-t), At(t-p)} provides a popula- 
tion of about 4.7X10 11 multi-valued sets and associated test 
25 subjects that can be discriminated. This number can be further 
increased by consideration of the cross-correlation values 
{(ap,aq), (ap,ar), (ap,as), (ap,at), (aq,ar), (aq,as), (aq,at), (ar, 
as), (ar,at), (as, at), (At(p-q),At(q-r)), (At(p-q),At(r-s)), (At(p- 
q),At(s-t)), (At(q-r),At(r-s)), (At(q-r),At(s-t)), (At(r-s),At(s-t)), 
30 (ap,At(p-q), (aq,At(p-q), (aq,At(q-r), (ar,At(q-r), (ar,At(r-s), 
(as,At(r-s), (as,At(s-t) or (at,At(s-t)} discussed in the preced- 
ing, so that the total number of test subjects that can be 
discriminated is estimated to be at least 1.2xl0 12 . 

Where the (normalized) mean value is used for comparison 
35 of candidate person parameters versus reference parameters, 
the ten statistically determined parameters in the set {ap, aq, 
ar, as, at, At(p-q), At(q-r), At(r-s), At(s-t), At(t-p)} can be 
compared in a visually perceptive display as illustrated in 
FIGS. 13 A and 13B. In the xy -graphs in FIG. 13 A, the mean 
40 peak amplitude values (y-axis, with the associated 
signums±), m„(ap;meas) vs. m„(ap;ref), m„(aq;meas) vs. 
m w (aq;ref), m w (ar;meas) vs. m w (ar;ref), m w (as;meas) vs. 
m„(as;ref), m„(at;meas) vs. m„(at;ref), and m„(ap;meas) vs. 
m„(ap;ref), are compared, in pairs, by vertical separation of 
45 horizontal line segments, one above the other for each pair. 
These pairs of peak amplitude values can be compared with 
reference to the respective y-axis values, m w (a;meas) and 
m w (a;ref) in FIG. 13 A. 

The x-axes, representing time, in FIG. 13A are aligned 
50 with the same initial value, t=t0(ref)=t0(meas), representing 
the times of the respective peak amplitudes for the P-wave, 
a=ap, and a final time value, t=tl(ref) or t=tl(meas), repre- 
senting the next succeeding time value for beginning of a new 
PQRST cycle for the reference complex; the value t=tl(ref) 
55 may not coincide with a corresponding time value, t=tl 
(meas) for the measured PQRST complex. 

A center point, corresponding to peak amplitude, for each 
horizontal segment in the sequence {m„(ap;meas), m„(aq; 
meas), m„(ar;meas). m„(as;meas), m„(at;meas), m„(ap; 
60 meas)} is separated in time (x-axis) from the next succeeding 
horizontal segment by a temporal distance m„(At;meas), with 
At=At(p-q), At(q-r), At(r-s), At(s-t) and At(t-p), in FIG. 13 A In 
a similar manner, a center point, corresponding to peak ampli- 
tude, for each horizontal segment in the sequence {m„(ap; 
65 ref), m„(aq;ref), m„(ar;ref). m„(as;ref), m„(at;ref), m„(ap; 
ref)} is separated in time (x-axis) from a center point for the 
next succeeding horizontal segment by a temporal distance 
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m„(At;ref). The pairs of pesk amplitude time intervals, m„(At 
(p-q) and m„(At(p-q);ref), m„(At(q-r) and m„(At(q-r);ref), 
m„(At(r-s) and m„(At(r-s);ref), m w (At(s-t) and m w (At(s-t) ;ref), 
and m w (At(t-p);meas) and m w (At(t-p);ref), can be compared 
with reference to the respective x-axis values, m w (At);meas) 5 
13 A and m w (At;ref) in FIG. 13 A). Display of these mean 
value pairs on the same graph may allow easier comparison of 
the corresponding mean values. The mean values, m w (a; 
meas) and m„(a;ref), can be replaced by the corresponding 
standard deviation values , sd w (a;meas) and sd w (a ;ref), in FIG. 

13 A, if desired. 

FIGS. 13C-13E graphically compares the mean values, 
m„(0;mode;meas) and m„(0;mode;ref), for the waves associ- 
ated with the atrial depolarization vectors V(0;ad;meas) and 
V(0;ad;ref), for the waves associated with the ventricular 
depolarization vectors V (0 ; vd;meas) and V(0 ; vd;ref), and for 1 5 
the waves associated with the ventricular repolarization vec- 
tors. V (0 ; vr;meas) and V (0 ; vr;ref), where each depolarization 
and repolarization vector V can be represented by an angle 0 
relative to an x-axis (0=0°) and a length L=l VI representing 
intensity of the vector V. The mean values, m w (0 ;mode;meas), 20 
m w (0;mode;ref), m w (L;mode;meas) and m w (L;mode;ref), can 
be replaced by the corresponding standard deviation values, 
sd w (0;mode;meas), sd M (0;mode;ref), sd M (L;mode;meas) and 
sd„(L;mode;ref), in FIGS. 13C-13E, if desired. In each 
instance, the “measured” value is a statistical value, measured 2 s 
over many heart cycles for the candidate person whose 
asserted identity is to be confirmed or refuted. 

Full application of the ECG analysis discussed here is, 
admittedly, fairly complex, although the richness and diver- 
sity of the results are often an ample reward for the effort. The 
analysis, results and conclusions discussed here can be 
applied wherever access to a facility or to selected informa- 
tion is limited to a small number of certified individuals. The 
applications include: access to information and/or to modifi- 
cation of information in banking and money handling/trans- 
fer activities; access to weapons storage and weapons main- 35 
tenance facilities; access to encryption codes and associated 
keys; and authentication (or refutation of) an identity asserted 
by a candidate person in a relatively short time interval (a few 
minutes). It may be possible, for some purposes, to replace 
this full ECG analysis with a modified version in which, as 40 
one example, timed measurements of signals sensed at two or 
more fingers or thumbs, pressed against separate sensing 
pads, and time shifted correlations are analyzed to distinguish 
between certain persons, activities; access to weapons storage 
and weapons maintenance facilities; access to encryption 45 
codes and associated keys; and authentication (or refutation 
of) an identity asserted by a candidate person in a relatively 
short time interval (a few minutes). It may be possible, for 
some purposes, to replace this full ECG analysis with a modi- 
fied version in which, as one example, timed measurements of 
signals sensed at two or more fingers or thumbs, pressed 50 
against separate sensing pads, and time shifted correlations 
are analyzed to distinguish between certain persons. 

Sensor placement, for purposes of reading signals that 
result in the PQRST complex may be constrained by religious 
or cultural values that limit or prohibit contact of any device 55 
with the human body. Where this occurs and comparison of 
the PQRST complexes is central to identity authentication, 
the candidate person may have to be denied access to what- 
ever site or object or information that person seeks. 

60 

APPENDIX 

Application of Maximum Likelihood Classification 
to Parameters 

65 

The asserted identity of a candidate person can be verified 
or refuted using a subset, or the full set, of the features dis- 
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cussed in the preceding, but employing a statistical classifi- 
cation method based on a Gaussian Mixture Model (GMM). 
The features, preferably in dimensionless format, can be con- 
catenated into a feature vector x of dimension D. Data from a 
particular ECG training set (reference person) are collected to 
construct a statistical model, referred to as a GMM. For 
verification purposes, new data are scored against the GMM 
that corresponds to the asserted identity. 

The verification uses a likelihood function, based on N 
mixture components, 

p(jt I Af)=2 f w f p f (jc), (A- 1 ) 

applied to a test feature vector x for the model M, where p.(x) 
is an ith mixture, modeled as a Gaussian density function, and 
w z is a non-negative weight value associated with the mixture 
number i and satisfying 

2^=1. (A-2) 

Each density function is assumed to have a Gaussian form, 

P.'( jc .'PiH( 23«) D,2 I|CJI‘' 2 }“ 1 exp "■(Cf 1 ) (x-fi,)/ 

2}, (A-3) 

where p. is a vector mean, determined from the training data 
and C z is a covariance matrix for the ith Gaussian mixture 
density for the model M, both determined from the training 
data (reference person), and x is a corresponding feature 
vector determined from the test data (candidate person). 
Training data are multiple feature vectors Xj recorded for the 
reference person at various times before testing of the candi- 
date person. 

The overall density function is assumed to be a product of 
the density functions p^(X^IM) corresponding to the jth vector 
density. A logarithm of the likelihood function 

L=\og{p{X\M)}=?.j log p(Xj\M) (A-4) 

is used here as a discriminat function, where X is a global 
vector referring to the collection of vectors {x 7 } 7 , and L is 
partially differentiated with respect to each of the vector 
parameters components p z to estimate optimum parameter 
values and any other parameters of interest. The solution 
parameters are substituted in the expression for p. in Eq. 
(A-3), together with measured values for the test vector x, and 
the density function numerical value (DFNV) is compared 
with a selected threshold value. If DFNV is at least equal to 
the threshold value, this condition is interpreted as indicating 
that the candidate person is likely the same as the reference 
person. If DFNV is less than the threshold value, this condi- 
tion may be interpreted as indicating that the candidate person 
is likely not the reference person. The MLE procedure is 
discussed in Jae Myung, “Tutorial on maximum likelihood 
estimation,” Jour. Of Math. Psychology, vol. 47 (2003), pp. 
90- 1 00, incorporated by reference herein. The MLE approach 
requires considerable effort to gather sufficient training data 
but has the advantage that use of an extraneous error metric, 
such as the Minkowski distance(s) in Eqs. (17)-(20), need not 
be imposed. 

Where optimal values of the vector means have been 
determined, using the procedure indicated above, the mea- 
sured vector values, or a mean value of the vector values x, 
can be inserted into Eq. (A-3), using the optimal values of the 
vector means p, z , and a numerical value for the density func- 
tion p(x;p.) can be computed. If this numerical value lies in a 
selected numerical range, this condition is interpreted as indi- 
cating that the andidate person is likely to be the reference 
person. If this numerical value does not lie in the selected 
numerical range, this condition is optionally interpreted as 
indicating that the andidate person is likely not the reference 
person. 
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What is claimed is: 

1. A method for estimating an identity of a candidate per- 
son, the method comprising: 

providing a sequence of measured values of at least one 
time interval length value, drawn from a set of measured 5 
time interval length values, At=At(p-q), At=At(q-r), 
At=At(r-s), At=At(s-t) and At=At(t-p), corresponding to a 
P wave, a Q wave, an R wave, an S wave and a T wave, 
respectively, for a sequence of PQRST electrical signal 
heart complexes for a candidate person; 
computing at least one measured statistical parameter, 
drawn from a set of parameters sp2(At;meas), compris- 
ing normalized mean m w (At;meas), normalized standard 
deviation sd w (At;meas), skewness sk(At;meas) and kur- 15 
to sis ku(At;meas), for the at least one measured time 
interval length, At=At(p-q), At=At(q-r), At=At(r-s), 
At=At(s-t) and At=At(t-p), for the sequence of PQRST 
complexes associated with the candidate person; 
providing at least one reference statistical parameter, 20 
drawn from a set sp2(At;ref), of reference parameters 
comprising normalized mean m w (At;ref), normalized 
standard deviation sd„(At;ref), skewness sk(At;ref) and 
kurtosis ku(At;ref), for at least one reference time inter- 
val length, At=At(p-q), At=At(q-r), At=At(r-s), At=At(s- 25 
t) and At=At(t-p), for PQRST complexes associated with 
a known reference person; 

forming a difference, A2(k2), numbered k2=l, . . . , K2 
(K2i^ 1), dependent upon a magnitude of at least one of 
the differences, m M ( At ;meas) -m„ (At ;ref), sd„(At; 30 

meas)-sd„(At;ref), sk(At;meas)-sk(At;ref), and ku(At; 
meas)-ku(At;ref), of statistical parameters drawn from 
the sets sp2(At;meas) and sp2(At;ref), multiplied by a 
non-negative weight value, wl', w2', w3* and w4’, 35 
respectively; 

when A2(k2) lies in a selected range R2(k2) of values, 
interpreting satisfaction of this condition as indicating 
that the candidate person is likely to be the reference 
person; and 40 

when A2(k2) does not lie in the range R2(k2) of values, 
interpreting non-satisfaction of this condition as indicat- 
ing that the candidate person is not likely to be the 
reference person; 

displaying or otherwise indicating at least one of (i) the 45 
interpretation that the candidate person is likely to be the 
reference person and (ii) the interpretation that the can- 
didate person is not likely to be the reference person. 

2. The method of claim 1, further comprising choosing said 

difference A2(k2) to comprise 50 

A2 (£2)={ wl' I m : (Ar;meas)-m„(Ar;ref) Z 2 + w2 ' \sd n (At; 
meas)-sd n (At;ref) [^+w3 ' ls&(Ar;meas)-s£(Ar;ref) 

\P+w4' \ku(At;mesLs)-ku(At;ref) Z 2 } llfl , 

where f2 is a selected positive number. 55 

3. The method of claim 1, further comprising: 
providing, as said reference person, a person who has an 

abnormal feature, associated with a medical condition 
that is not transitory, and associated with at least one of 
said at least one reference statistical parameters, for at 60 
least one of said time interval length values, At=At(p-q), 
At=At(q-r), At=At(r-s), At=At(s-t) and At=At(t-p); and 
when said candidate person has a measured statistical 
parameter, for said corresponding time interval length 
value At that does not include the abnormal feature, 65 
interpreting this condition as indicating that said candi- 
date person is likely not said reference person. 
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4 . The method of claim 1 , further comprising: 

providing, as said reference person, a person who does not 

have an abnormal feature, associated with a medical 
condition that is not transitory, associated with at least 
one of said at least one reference statistical parameters, 
for at least one of said time interval length values, At=At 
(p-q), At=At(q-r), At=At(r-s), At=At(s-t) and At=At(t-p), 
and 

when said candidate person has a measured statistical 
parameter, for said corresponding time interval length 
value At that does include the abnormal feature, inter- 
preting this condition as indicating that said candidate 
person is likely not said reference person. 

5. The method of claim 1, further comprising: 

displaying, at least one of said measured statistical param- 
eter values, normalized mean m M (At;meas), normalized 
standard deviation sd„(At;meas), skewness sk(At;meas) 
and kurtosis ku(At;meas); and 

displaying at least one of said reference statistical param- 
eter values, normalized mean m„(At;ref), normalized 
standard deviation sd w (At;ref), skewness sk(At;ref) and 
kurtosis ku(At;ref) corresponding to the at least one of 
said measured statistical parameter values. 

6. The method of claim 5, further comprising displaying 
said at least one of said measured statistical parameter values, 
and said least one of said corresponding reference parameter 
values, on a first xy-graph. 

7. The method of claim 6, wherein: 

said at least one measured statistical parameter value is 
displayed on said first graph as a line segment, substan- 
tially parallel to an x-axis of said first graph with a first 
y- value on a y-axis of said first graph corresponding to 
said at least one measured statistical parameter value; 
and 

said at least one corresponding reference statistical param- 
eter value is displayed on said first graph as a line seg- 
ment substantially parallel to the x-axis on said graph, 
with a second y-value on the y-axis of said first graph 
corresponding to said at least one corresponding refer- 
ence statistical parameter value, where the x-axis of said 
first graph represents time for said measured statistical 
parameter values, At=At(p-q), At=At(q-r), At=At(r-s), 
At=At(s-t) and At=At(t-p), and represents time for said 
reference statistical parameter values, At=At(p-q), At=At 
(q-r), At=At(r-s), At=At(s-t) and At=At(t-p). 

8. The method of claim 5, further comprising displaying 
said at least one measured statistical parameter value, on a 
first xy-graph, and displaying said at least one corresponding 
reference parameter value on a second xy-graph. 

9. The method of claim 6, further comprising: 

displaying said at least one measured statistical parameter 

value on said first graph as a line segment substantially 
parallel to an x-axis on said first graph, with a first 
y-value on a y-axis of said first graph corresponding to 
said at least one measured statistical parameter value; 
and 

displaying said at least one corresponding reference statis- 
tical parameter value on a second graph as a line segment 
substantially parallel to an x-axis on said second graph, 
with a y-value on a y-axis of the second graph corre- 
sponding to said at least one reference statistical param- 
eter value, where an x-axis of said first graph, represent- 
ing time, and an x-axis of the second graph, representing 
time, are aligned with reference to (i) a time of occur- 
rence of said measured time interval length value, At=At 
(p-q), At=At(q-r), At=At(r-s), At=At(s-t) and At=At(t-p), 
for said measured statistical parameter values and (ii) a 
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time of occurrence of said corresponding reference peak 
amplitude value, At=At(p-q), At=At(q-r), At=At(r-s), 
At=At(s-t) and At=At(t-p), for said reference statistical 
parameter values. 



